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of 471 x 10 6 BB pairs collected by the BABAR experiment. We fully reconstruct the hadronic decay 
of one of the B mesons in the T(4S) — > BB decay, and search for the B —>■ K^*'vv decay in the 
rest of the event. We observe no significant excess of signal decays over background and report 
branching fraction upper limits of B{B+ ~¥ K+vv) < 3.7 x lCT 5 , B(B° ->■ K°vv) < 8.1 x 1(T 5 , 
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INTRODUCTION 

Flavor-changing neutral-current (FCNC) transitions, 
such as b — > svv, are prohibited in the Standard Model 
(SM) at tree-level. However, they can occur via one-loop 
box or electroweak penguin diagrams, as shown in Fig.[T] 
They can occur also in the SM via a quarkonium reso- 
nance state b — > sec, cc — > vv, where the cc decay is me- 
diated by a virtual Z° boson (Fig. [2]) . This latter decay 
process has the same final state as b — > svv with an addi- 
tional constraint from the on-shell cc mass. Both the b — >• 
svv and cc — > vv decay rates are expected to be small 
within the SM, with branching fractions estimated to be 
B{B+ -> K+vv) = B(B° -> K°vv) = (4.5 ± 0.7) x 1(T 6 , 
B(B+ K*+vv) = B(B° -> K*°vV) = (6.8±^;?) x lO" 6 
|, and B(J/tp -> vv) = (4.54 x 10- 7 ) • B(J/ip' -> e+e~) 
3|. The b — > svv rates are predicted with smaller the- 
oretical uncertainties than those in the corresponding 
b — > s£ + £~ modes due to the absence of long-distance 
hadronic effects from electromagnetic penguin contribu- 
tions. 

Various new-physics scenarios exist that could signif- 
icantly enhance the b — > svv branching fractions, as 
well as modify the expected SM decay distributions of 
sb = q 2 /m B , where q 2 is the squared magnitude of 
the four-momentum transferred from the B meson to 




u, c, t 



u, c, t 



FIG. 1: Lowest-order SM Feynman diagrams for b — > svv 
transitions. The virtual top quark provides the dominant con- 
tribution in each case. 



the neutrino pair, and is the B meson mass. Some 
of these scenarios predict massive particles that could 
contribute additional loop diagrams with similar ampli- 
tudes as those in the SM, such as nonstandard Z° cou- 
plings with supersymmetric (SUSY) particles [l|, fourth- 
generation quarks Q, anomalous top-charm transitions 
\m, or a massive U(l) gauge boson Z' 0, S- bmce 
b — > svv has two final-state neutrinos, other sources of 
new physics can also contribute to the experimental sig- 
nature of a kaon and missing four-momentum, such as 
low-mass dark- matter (LDM) candidates 0, > unpar- 
ticles @, right-handed neutrinos @, or SUSY particles 
[icj | . Models with a single universal extra dimension also 



11] 



predict higher decay rates 

The decays Jjif) — > vv and ip(2S) — > vv pro- 
vide additional windows for new-physics searches. In 
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FIG. 2: Lowest-order Feynman diagrams of (from left to 
right) the SM decay cc — > vv, the SUSY decay cc into a pair of 
goldstinos (<?) via a c-squark in the t-channel, and the SUSY 
decay cc — > gg via a virtual Z° in the s-channel. 



spontaneously-broken SUSY, a cc resonance can decay 
into a pair of goldstinos via either a virtual Z° in the 
s-channel or a c-squark exchange in the t-channel [2} 
(Fig. [2J . The contribution of a massive SU(2) gauge bo- 
son Z', introduced in the left-right SUSY model, could 
suppress the decay rates up to an order of magnitude [1| ■ 
Conversely, a low-mass U(l) gauge boson U could en- 
hance the invisible decay rates of quarkonium states by 
several orders of magnitude by coupling to LDM particles 
[HI, 0|. The U boson could decay into a pair of spin- 1 /^ 
Majorana (xx)> spin- 1 /^ Dirac (xx)> 01 spin-0 (<p<p) LDM 
particles. 

We search for B — > KvV and B — > K*vv, and for 
J/V» — !• z^ 7 and ^(25) -> vv via 5 -> K^J/i/j and 
£ — > K(*)ip(2S) respectively, where signifies a 

charged or neutral K or X* meson [3]. We use a tech- 
nique in which one B meson is exclusively reconstructed 
in a hadronic final state before looking for a signal decay 
within the rest of the event. Since the four-momentum of 
one B meson is fully determined, the missing mass reso- 
lution on the two final-state neutrinos and the suppres- 
sion of background are improved with respect to other 
reconstruction techniques. 

Several previous searches for B — > KvV and B — > 
K*vv have been performed by both the BABAR and 
BELLE collaborations 

Currently, the most str n- 
gent published upper limits at 90% confidence level (CL) 
are B(B+ -> K+vv) < 1.3 x lfr 5 [HI and B(B -> 
ifVF) < 8 x 10~ 5 18]. The -> limit was 

determined using semileptonic-tag reconstruction, which 
produces samples that are statistically larger and inde- 
pendent of those produced using the hadronic-tag recon- 
struction employed in this search. The B(B —> K*vv) 
limit was a combination of two BABAR analyses, one us- 
ing semileptonic-tag reconstruction and the other using 
hadronic-tag reconstruction. 

A J/tp — > vv search via ip(2S) — > rnrJ/tp was performed 
by the BES collaboration, which set an upper limit at 
90% CL of B(J/ip -> vv) < 1.2 x 10~ 2 • B(J/ip -> fi + fi~) 
fljj ]. This article presents the first search for J/ip — > 
vv using the hadronic-tag reconstruction of a B meson 
decay. A search for ip(2S) — > vv has not been performed 
previously. 



THE BABAR DETECTOR AND DATA SAMPLE 



This search uses a data sample of 471 ± 3 million 
BB pairs, corresponding to an integrated luminosity of 
429 fb" 1 collected at the T(AS) resonanceEJ. The data 
were recorded with the BABAR detector [21( at the PEP- 
II asymmetric-energy e + e _ storage rings. The charged- 
particle tracking system consists of a five-layer double- 
sided silicon vertex tracker and a 40-layer drift cham- 
ber, both coaxial with a 1.5 T solenoidal magnetic field. 
Charged kaons and pions are distinguished by specific 
ionization energy-loss measurements from the tracking 
system for lower momentum particles, and by measure- 
ments from a ring-imaging Cherenkov radiation detector 
for higher momentum particles. A CsI(Tl) electromag- 
netic calorimeter is used to reconstruct photons of energy 
greater than 20MeV and to identify electrons. Muon 
identification is provided by the instrumented flux re- 
turn of the magnet. Particle identification (PID) algo- 
rithms are trained to identify charged particle types by 
using 36 input parameters including momentum, polar 
and azimuthal angles, the Cherenkov angle, and energy- 
loss measurements [22j. We employ PID criteria that 
select K + mesons with an efficiency greater than 85% 
and with approximately 1% misidentification probability 
for pions and muons. 



Signal and background decays are studied using Monte 
Carlo (MC) samples simulated with Geant4 (23j. The 
simulation includes a detailed model of the BABAR detec- 
tor geometry and response. Beam-related background 
and detector noise are extracted from data and are 
overlaid on the MC simulated events. Large MC sam- 
ples of generic BB and continuum (e + e~ — > t + t~ or 
e + e~ — > qq, where q — u,d, s, c) events provide ten times 
the number of T(4S) — > BB and e + e~ — > cc events as in 
the data sample, and four times the number of other con- 
tinuum decays. The T(45) — > BB signal MC samples are 
generated with one B meson decaying via vv, 
with and without the cc resonances, while the other B 
meson decays according to a model tuned to world av- 
erages of allowed decay channels. The sb distributions 
for B —> K^*'vv decays within signal MC samples are 
generated initially using a phase-space model, and then 
reweighted using the model from Ref. [l|, henceforth re- 
ferred to as ABSW. Within B -» K*vv decays, this 
model is also used to reweight the helicity-angle distri- 
bution between the signal B and the K + or K° flight 
directions in the K* rest frame. The helicity amplitudes 
for the decay channels B — > K*J/ij) and B — > K*ip{2S) 
are generated using values taken from a BABAR measure- 
ment [24 1 . 
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ANALYSIS METHOD 

Event selection for both the B — >• K^vV and B -» 
cc, cc — > vv searches begins by fully reconstructing 
a B meson (Stag) in one of many hadronic final states, 
B -> SX had , where S is a "seed" meson (£>«+, flW , 

£>i* )+ , or J/V>) and Xr, is a collection of at most five 
mesons, composed of charged and neutral kaons and pi- 
ous with a net charge of —1. This method, which was 
used also in Ref . [25[ , reconstructs additional modes with 
resp ect to previous hadronic-tag B —> K^yy analyses 
jl5l [l8j ]. and results in approximately twice the recon- 
struction efficiency. The D seeds are reconstructed in 
the decay modes D + K°n + , K°tt^ 



x10 J 



K~ 



K' 



' 7T 

K~ 



K-jt 

K'TT 

K + K~, 7T + 7T~, 7r + 7T~7r° 

are reconstructed as D* + 
D° r , D*+ -> £+ 7 ; Dr 
and J/V> — > e + e , 
constructed via their decay to 7r + 7r~ 



'°7r+7r°, ^°7r+7r+7r- 
if+k-7r+, K+K-ir+ir : D° 



and K°ir°. 



— ► 

"7T - , i^OTT+TT-TT , 

Additional seeds 
vr ; D*° -> Z^ , 



The if 9 candidates are re- 



Well-reconstructed S 



two kinematic variables: A.E = Eb^ 
where i?B ta , and p Bta 



tag candidates are selected using 
— y/s 1 2 and Toes = 
are the energy and 



momentum vector of the B tag candidate, respectively, in 
the e + e~ center-of-mass (CM) frame and y/a is the total 
energy of the e + e~ system. The value of AE, which 
peaks at zero for correctly reconstructed B mesons, is 
required to be between —0.12 and 0.12 GeV or within two 
standard deviations around the mean for a given A^~ ad 
mode, whichever is the tighter constraint. If more than 
one St a g candidate is reconstructed, the one in the mode 
with the highest purity (fraction of candidates that are 
correctly reconstructed within a given B tag decay mode) 
is chosen. If there are multiple candidates with the same 
purity, the one with the smallest \AE\ is selected. 

After requiring that the event contains between one 
and three tracks not used in the B tag reconstruction 
("signal-side" tracks), the purity of each mode is recalcu- 
lated, and only the B tag modes that have a recalculated 
purity greater than 68% are retained. This results in a 
total of 448 final states. This purit y va lue was optimized 
by maximizing the figure of merit [26j 



(1) 



N, 



bko 



where the number of sigmas n a — 1.28 corresponds to a 
one-sided Gaussian limit at 90% CL, ef 6 is the total sig- 
nal efficiency, and N^ ks is the expected number of back- 
ground events, with i representing one of the signal decay 
channels. All other selection criteria discussed henceforth 
were optimized simultaneously using this same figure of 
merit. 
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FIG. 3: The tties distribution for the Stag candidates in data 
(points) and in the expected combinatorial background as pre- 
dicted by the MC (shaded). This distribution includes only 
the charged and neutral Stag candidates that pass the pu- 
rity restrictions, the multivariate continuum suppression, and 
a requirement of one to three signal-side tracks. The data 
within the J71es sideband region is used to extrapolate the 
expected number of combinatorial background events within 
the signal region. 



The signal region of the -Btag candidate is defined as 
5.273 < m ES < 5.290 GeV/ c 2 (Fig. [3]), since correctly re- 
constructed B mesons produce a peak in this region near 
the nominal _B-meson mass. The -Btag candidates that 
are incorrectly reconstructed ("combinatorial" events), 
which result from continuum events or are due to tracks 
or clusters assigned to the wrong B meson, produce a dis- 
tribution that is relatively uniform below the toes signal 
region and decreases toward the kinematic limit within 
it. Approximately 0.3% of signal MC events and 12.0 
million data events contain a -B tag that is reconstructed 
using the above requirements and found to be within the 
Toes signal region. 

Since B mesons are spin zero and are produced with 
low momentum in the CM frame (~ 0.32GeV/c), their 
decay products are more isotropically distributed than 
non- .B-B background. For example, |cos0t|, where 0j is 
the angle in the CM frame between the B tag thrust [27J 
axis and the thrust axis of all other particles in the event, 
has a uniform distribution for BB events but peaks near 
one for continuum events. Continuum background is sup- 
pressed by using a multivariate likelihood selector based 
on six event-shape variables. These consist of |cos 9t\ , the 
cosine of the angle between p*B tag and the beam axis, the 
magnitude of the -B tag thrust, the component of the -B tag 
thrust along the beam axis, the angle between the miss- 
ing momentum vector (p m i ss ) and the beam axis, and the 
ratio of the second-to-zeroth Fox- Wolfram moment [28| 
computed using all charged and neutral particles in the 
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event. The multivariate selector requires 

where V q (xj) and Vsixj) are probability density func- 
tions determined from MC that describe continuum and 
signal-like BB events, respectively, for the six event- 
shape variables Xj. The Lb requirement, which was op- 
timized with other selection criteria using Eq. ([1]), also 
improves the agreement between data and MC by sup- 
pressing unmodeled continuum backgrounds. 

In the sample of selected B tag candidates, signal events 
are chosen such that all of the remaining tracks can 
be combined to yield a sing le K& candidate. The 
sum of the K^*' and B tag candidate charges must equal 
zero. Since signal decays have two final-state neutri- 
nos, these events are required to have missing energy 
greater than zero, where the missing energy is defined 
as the CM energy minus all detected calorimeter de- 
posits from charged and neutral particles in the event. 
For B — > K(*>vV, the signal decays are reconstructed 
in six channels: B + — > K + vV\ B° —> K a vv where 
K° -> K° s - B+ -> K* + uV, where K*+ -> K+ir° and 
K*+ -> K°ir+; and B° -> K*°vV, where K*° -> K+tt~ 
and K*° —> K^tt . For cc — > i/V, the same six signal 
channels are employed with an additional requirement 
that the momentum is consistent with a two-body 
decay, cither B -> K^J/ip or B -> K^ip(2S). The J/ip 
and tp(2S) mesons then decay into a pair of neutrinos, 
thus yielding the same final states as for B — > K^*'vV. 

We reconstruct K® — > n + n~ decay candidates using 
two tracks of opposite charge, which originate from a 
common vertex and produce an invariant mass within 
±7MeV/c 2 of the nominal K° s mass 29]. The PID for 
each track must be inconsistent with that for an electron, 
muon, or kaon. The ir° candidates are reconstructed 
from pairs of photon candidates with individual ener- 
gies greater than 30MeV, a total CM energy greater 
than 200 MeV, and a 77 invariant mass between 100 and 
160MeV/c 2 . All K+ candidates must satisfy the PID 
criteria for a kaon. 

Reconstructed K* candidates are required to have an 
invariant mass within ±70 MeV/ c 2 of the nominal K* 
mass 29]. A K* + — > K®tt + candidate combines a Kg 
candidate with a track that satisfies the PID criteria for 
a pion. If more than one K*+ -> K°n+ candidate can 
be reconstructed in an event, the one with the mass clos- 
est to the nominal K* + mass is chosen. A K*° — > K + ir~ 
candidate combines one track that satisfies the PID crite- 
ria for a kaon with one that is inconsistent with the PID 
criteria for an electron, muon, or kaon. In an event con- 
taining a K + (Kg) candidate and no additional signal- 
side tracks, K* — > K + tt° (K^tt ) candidates are recon- 
structed if the invariant mass of a 7r° candidate and the 
K + (Kg) candidate falls within the K* mass window; 



otherwise the event is considered for the K + (Kg) signal 
channel. If more than one K*+ -> K+ir° or K*° -> K°tt° 
candidate can be reconstructed, the one with the highest 
energy n° candidate is chosen. 

Once the B ta , g and K^ are identified, the signal events 
are expected to contain little or no additional energy 
within the calorimeter. However, additional energy de- 
posits can result from beam-related photons, hadronic 
shower fragments that were not reconstructed into the 
primary particle deposit, and photons from unrecon- 
structed D* — > Dj/tt° transitions in the B tag candidate. 
Only deposits with energy greater than 50 MeV in the 
rest frame of the detector are considered, and the sum 
of all such additional energy deposits (E ex t ia ) is required 
to be less than a threshold value (Ei). The values of Ei, 
given in Table [H were optimized with the other selection 
criteria but were allowed to differ between signal chan- 
nels. For events within the K + signal channel, clusters 
identified as kaon shower fragments are not included in 
the -Eextra sum. A fragment candidate is defined as a 
cluster whose momentum vector, when compared to that 
of the signal track, is separated by polar and azimuthal 
angles (relative to the beam axis and in the rest frame 
of the detector) of A8 and Acf>, respectively, such that 

r clus < 15°, where r clus = (A9) 2 + |(Q K • A0 - 8°) 2 

and Qk = ±1 is the K^ charge. The r c i us and frag- 
ment candidate definitions were optimized using studies 
of truth information in the signal MC samples. The re- 
covery of these kaon shower fragments improves the final 
signal efficiency in the K + channel by about 13%. This 
procedure was explored for the other signal tracks, but 
the effect was small. 



TABLE I: Threshold values Ei for the EWtra variable in each 
of the signal channels, determined using Eq. |l]). The channels 
in brackets refer to the K* decay products. 



Channel K+ 


K° 


[K+ir ] 


[K s ir + ] 


[K+TT-] 


[^°] 


Ei [GeV] 0.11 


0.28 


0.18 


0.29 


0.31 


0.33 


The searches 


for B 




vis and 


for cc - 


> vv via 



B — > K^cc diverge in the final step of the signal selec- 
tion, which involves restricting the kinematics of the de- 
cay. The value of sb is calculated as (ps 3ig —Pkm ) 2 /m B , 
where p^M is the four-momentum of the K^ candidate, 
and PB Bis is the expected signal B four-momentum with 
an energy of \fs/2, the nominal _B-meson mass, and a 
momentum vector pointing opposite the Stag momen- 
tum. For B — > K^vV, the signal region optimized for 
maximum SM sensitivity is < sb < 0.3 for all six signal 
channels. This corresponds to a K^ momentum greater 
than about 1.8 (1.7) GcV/c in the signal B rest frame 
for B —> KvV (B — > K*vv) events. Partial branch- 
ing fractions over the full sg range are also provided 
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for sensitivity to new-physics scenarios that modify the 
kinematic distributions for B — > K^*'vv, For cc — > vv 
via B — > K^cc, the invariant mass of the two neutri- 
nos m vV = y/ ' SBin 2 B is expected to correspond to the 
mass of the J/tp (3.097 GeV/c 2 ) meson or to that of the 
ip(2S) (3.686 GeV/c 2 ) meson. Signal events are selected 
within three standard deviations around the nominal cc 
masses, which results in windows of 3.044 < muu < 3.146 
(3.019 < m vV < 3.175) GeV/c 2 for the B -> KJ/ip 
(B ->■ K*J/i>) channel, and 3.650 < m vV < 3.724 
(3.627 < m vV < 3.739) GeV/c 2 for the B -> Kip(2S) 
(B ->■ K*ip(2S)) channel. 

To avoid experimenter bias, all the above selection cri- 
teria and values were optimized using the MC before 
looking at any data events within the -E extra and Toes 
signal regions. 



BACKGROUND AND BRANCHING FRACTION 
EXTRACTION 

The total number of background events s in the 
signal region has two components: iV I peak is the num- 
ber of expected background events having a correctly 
reconstructed Stag candidate and hence peaking within 
the Toes signal region, and Nf omh is the number of ex- 
pected combinatorial background events, including both 
continuum events and BB events with an incorrectly re- 
constructed -Btag candidate. To reduce the dependence 
on MC simulations, the number of Nf olnh events is ex- 
trapolated directly from the observed data events within 
the toes sideband region, defined as 5.200 < tties < 
5.265 GeV/c 2 and depicted in Fig. |3l The shape of the 
combinatorial toes distribution is estimated using MC 
samples of continuum events and of BB events recon- 
structed with the wrong charge. 

The number of Nf cak events is estimated from generic 
BB MC samples. Over half of Nf eak is found to be from 
B — » D(*'£v (I = e or /i) decays in which no lepton 
candidate is identified in the event and the K^*> is a 
daughter of the D or D* meson. One particular peaking 
background in the B — > K^vv search is B + t + v t , 
with t + — > K^ + v T , which has the same final state as the 
signal decay [3(| . Exclusive B + — > t + v t MC samples, 
assuming a branching fraction of (1.65 ± 0.34) x 10~ 4 
[29j . indicate that this background constitutes less than 
15 (5)% of the total background in the B + — > K + vv 
(B + K* + vV) channel. 

Since both Nf cak and ef s are determined from MC 
samples, we normalize the MC yields to the data to ac- 
count for differences between data and MC, such as from 
the Stag reconstruction and the modeled branching frac- 
tions of Stag modes within the MC. This normalization 
is performed before applying the full signal selection in 
order to have a large background-to-signal ratio; looser 
mass windows and -E C xtra selection requirements are 



used such that the number of background events is ap- 
proximately 60 times larger than the final background 
contribution, over the full sg spectrum. The peaking 
background component in the BB MC is then normal- 
ized to the number of data events that peak within the 
mES signal region. This peaking yield normalization is 
performed separately for charged and neutral -Btag can- 
didates, and results in the scaling of all signal and back- 
ground MC samples by 1.027 ± 0.039 (1.017 ± 0.044) for 
charged (neutral) -Btag candidates. 

The signal branching fractions are calculated using 



Bi 



N 



obs 



peak 



n: 



comb\ 



N 



(3) 



BB 



where N BB - = 471 x 10 6 is the total number of B me- 
son pairs in the data sample and N° hs is the number of 
data events within the signal region. The total signal ef- 
ficiency includes that of the B tag reconstruction and 
is determined separately for each of the signal channels 
i. Since misreconstructed events from other signal chan- 
nels contribute to Nf cak , the branching fractions of all 
signal channels are determined simultaneously. Branch- 
ing fraction limits and uncertainties are computed us- 
ing a mixed frequentist-Bayesian approach described in 
Ref . [3l| , with the systematic uncertainties on N^ ks and 
e^ lg modeled using Gaussian distributions. To combine 
the results of signal decay channels, we find the Bi value 
that maximizes a likelihood function defined as the prod- 
uct of the Poisson probabilities of observing iV° bs events. 



SYSTEMATIC STUDIES 



To verify the modeling of ef s and N i s , a control 
sample of B — > Dlv events is selected. In place of 
a signal K* candidate, the events are required to con- 



tain a reconstructed D° — > K 7r + , D~ 



K+TT- 



candidate with an invariant mass within 



±35 McV/c 2 of the nominal D -meson mass values [291 . 
The event must have one additional track that satis- 
fies the PID criteria of either an electron or muon. All 
other reconstruction and signal selection requirements 
are retained. The resulting yields in the data agree with 
MC expectations, assuming the well-measured branching 
fractions of B ->• Dlv within the 7% (12%) statisti- 
cal uncertainty of the data in the < Sb < 0.3 (J/ip or 
Tp{2S) mass) region. 

The control sample is used to determine the system- 
atic uncertainties due to the MC modeling of the -E ex tra 
variable within data. Additional uncertainties on jVP eak 
and are due to the K® and K* mass reconstruction 
windows, the ir° reconstruction, and the uncertainties in 
the branching fractions [29} of the dominant backgrounds 
contributing to Nf eak . The uncertainty on Nf omh is dom- 
inated by the sideband data statistics. Other systematic 
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uncertainties, such as those from PID, tracking, -Btag re- 
construction, N B -g, and the assumption that charged and 
neutral BB pairs are produced at equal rates, are all 
accounted for by the normalization of the MC peaking 
yields. Because the peaking yield in data depends on 
the extrapolated shape of the combinatorial B tag back- 
ground, the normalization scale factors are re-evaluated 
by varying the method used to extrapolate this shape. 
The resulting variations on the final N^ ks and s val- 
ues are taken as the systematic uncertainies due to the 
normalization. 

Due to the approximately 1.0% resolution on the sb 
measurement around sb — 0.3, an uncertainty is eval- 
uated within the B — > K^*'W signal region. Similarly, 
the resolution on m v v contributes to uncertainties within 
the J/ip — > vv and ip(2S) — > vv signal regions. Only the 
systematic uncertainties due to the iV 4 pcak branching frac- 
tions and to sb or m v -v differ between the B — > K^vv, 
J/ip — > vv, and ip(2S) — > vv searches. The systematic 
uncertainties are summarized in Tables [TTI and IIII1 the for- 
mer lists the uncertainties shared by the searches, while 
the latter lists those that differ. 



RESULTS FOR B 



Figure [4] shows the observed data yields, expected 
background contributions, and SM signal distributions 
over the full sb spectrum. Tables IIVI and [V] summa- 
rize the number of observed data events within the sb 
signal region (0 < sb < 0.3), expected backgrounds, 
B K^vv signal efficiencies, branching fraction cen- 
tral values, and branching fraction limits at the 90% CL. 
Combining the signal channels, we determine upper lim- 
its of B(B ->■ Kvv) < 3.2 x 10~ 5 and B{B K*vV) < 
7.9 x 10 -5 . Since we see a small excess over the expected 
background in the K + channel, we report a two-sided 
90% confidence interval. However, the probability of ob- 
serving such an excess within the signal region, given 
the uncertainty on the background, is 8.4% which cor- 
responds to a one-sided Gaussian significance of about 
1.4cr. Therefore, this excess is not considered significant. 

Using the same procedure as when combining signal 
decay channels, the B — ¥ Kvv branching fraction cen- 
tral values are combined with a previous semileptonic-tag 
BABAR analysis that searched within a statistically inde- 
pendent data sample 17J. We obtain combined EABAR 
upper limits at the 90% CL of 



B{B+ -> K + vv) < 1.6 x 10~ 5 , 
B(B° K°vv) < 4.9 x 10~ 5 , 
B(B -» Kvv) < 1.7 x 10~ 5 . 



and 



(4) 



The combined central value is B(B — » Kvv) = 
(0.8^0 6) x 10 -5 , where the uncertainty includes both 
statistical and systematic uncertainties. These combined 



TABLE II: Summary of systematic uncertainties that are 
shared by the B —¥ K^vv, J/tp — > vv, and ip(2S) — > vv 
searches. All values are relative uncertainties in %. The chan- 
nels in brackets refer to the K* decay products. 



Source 



K+ [K + TT°] [K°TT+] K° [K+TT-] [K°TT° 



ej lg normalization 3.5 3.5 

iV* kg normalization 2.3 2.3 
Kg reconstruction 

K* reconstruction 2.8 

7T° reconstruction 3.0 

£cxtra 4.5 6.0 



3.5 
2.3 
1.4 
2.8 



8.9 
6.0 
1.4 



8.9 
6.0 

2.8 



6.5 6.0 6.0 



6.0 
1.4 
2.8 
3.0 
6.5 



TABLE III: Summary of systematic uncertainties that dif- 
fer between the B — > K^W, J/ift — s> z/i7, and ip(2S) — s> z/i7 
searches, and the total systematic uncertainties for each sig- 
nal channel. All values are relative uncertainties in %. The 
total systematic uncertainties are determined by adding in 
quadrature each relevant uncertainty, including those listed 
in Table [Til 



Source 


K+ 


[K+n°] 




K° 


[K+TT-] 


[K° s n°] 






B -> K 










Aff eak B's 


2.8 


2.8 


2.8 


2.8 


2.8 


2.8 


sb resolution 


3.6 


3.6 


3.6 


3.6 


3.6 


3.6 


Total Nf e ^ syst. 


6.8 


8.9 


8.8 


9.7 


10.0 


10.9 


Total AT mb syst. 


2.3 


2.3 


2.3 


6.0 


6.0 


6.0 


Total ef g syst. 


6.7 


8.8 


8.8 


11.4 


11.7 


12.4 








vv 








jypeak g , g 


3.5 


3.5 


3.5 


3.5 


3.5 


3.5 


m„- resolution 


1.1 


2.1 


0.4 


0.7 


0.3 


1.3 


Total 7Vf eak syst. 


6.2 


8.6 


8.4 


9.3 


9.6 


10.5 


Total AT mb syst. 


2.3 


2.3 


2.3 


6.0 


6.0 


6.0 


Total ef g syst. 


5.8 


8.3 


8.0 


10.8 


11.1 


11.9 






V(2S) - 


-> vv 








Nf eak B's 


2.8 


2.8 


2.8 


2.8 


2.8 


2.8 


rrw resolution 


0.8 


2.4 


1.0 


0.9 


1.8 


3.1 


Total Nf e ^ syst. 


5.8 


8.5 


8.1 


9.1 


9.5 


10.7 


Total AT mb syst. 


2.3 


2.3 


2.3 


6.0 


6.0 


6.0 


Total ef s syst. 


5.8 


8.4 


8.1 


10.9 


11.2 


12.2 



results reweight the sb distribution to that of the ABSW 
theoretical model (dashed curve in Fig. |4|), which de- 
creases the signal efficiencies published in Ref. pjj by 
approximately 10%. The B — > K*vv central values also 
can be combined with the semileptonic-tag results from 
a previous BABAR search [18] . In order to obtain approxi- 
mate frequentist intervals, the likelihood functions in the 
previous search are extended to include possibly negative 
signals. We obtain combined BABAR upper limits at the 
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TABLE IV: Expected B — > K*vV background yields iV* kg = iVf cak + N£ omb , signal efficiencies ef s , number of observed data 
events iV° bs , resulting branching fraction upper limits at 90% CL, and the combined upper limits and central values, all within 
the < sb < 0.3 region. Uncertainties are statistical and systematic, respectively. The channels in brackets refer to the K* 
decay products. 





B+ -t [K+7r°]vV 


B+ -> [K° s ir + ]vV 


5° -> [K+tt~]vV 


B° -> [K° s ir ]vV 


^ypcak 


1.2 ± 0.4 ± 0.1 


1.3 ± 0.4 ± 0.1 


5.0 ± 0.8 ± 0.5 


0.2 ± 0.2 ± 0.0 


^ycomb 


1.1 ± 0.4 ± 0.0 


0.8 ± 0.3 ± 0.0 


2.0 ± 0.5 ± 0.1 


0.5 ± 0.3 ± 0.0 


AT bkg 


2.3 ± 0.5 ± 0.1 


2.0 ± 0.5 ± 0.1 


7.0 ± 0.9 ± 0.5 


0.7 ± 0.3 ± 0.0 


ef g (xlO- 5 ) 


4.9 ± 0.2 ± 0.4 


6.0 ± 0.2 ± 0.5 


12.2 ± 0.3 ± 1.4 


1.2 ± 0.1 ± 0.1 


N obs 


3 


3 


7 


2 


Limit 


< 19.4 x 10 -5 


< 17.0 x 10" 5 


< 8.9 x 10~ 5 


< 86 x 10~ 5 




fo o+6.2+1.7\ 
\°- -3.6-1.3) 


x 10~ s 


(2.Q±l 2 3 ±i° 7 ) x lO" 5 


Limit 


< 11.6 x 10 -5 


< 9.3 x 10" 


-5 


B(B ->■ K*vV) 




(2.7^:2 


-to) x io- 5 




Limit 




< 7. 


3 x 10~ 5 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 4: The sb distribution for (from top to bottom) B + — > 
K+W, B° -> K°vV, B+ ->■ K* + uV, and 5° -> K*°vv events 
after applying the full signal selection. The expected combi- 
natorial (shaded) plus mss-peaking (solid) background con- 
tributions are overlaid on the data (points). The signal MC 
distribution (dashed) is normalized to a branching fraction 
of 20 x 10~ 5 for B + -> K + vV and 50 x 10~ 5 for the other 
channels. Events to the left of the vertical lines are selected 
to obtain SM-sensitive limits, while the full spectra are used 
to determine partial branching fractions. 



TABLE V: Expected B -> KvV background yields N? ks = 
jypea. _j_ jy? omb ( signal efficiencies e° lg , number of observed 
data events N° bB , resulting branching fraction upper limits 
at 90% CL, the central values Bi, and the combined upper 
limits and central value, all within the < sb < 0.3 region. 
Lower limits at 90% CL are also reported, as discussed in 
the text. Uncertainties are statistical and systematic, respec- 
tively. The B° -> K°W efficiency accounts for B{K° -> K%) 
and B(K% -> tt+tt") @. 





B+ -> K+vV 


B° -> K°vV 


jypcak 


1.8 ± 0.4 ± 0.1 


2.0 ± 0.5 ± 0.2 


jycomb 


1.1 ± 0.4 ± 0.0 


0.9 ± 0.4 ± 0.1 




2.9 ± 0.6 ± 0.1 


2.9 ± 0.6 ± 0.2 


4 s (xlO" 5 ) 


43.8 ± 0.7 ± 3.0 


10.3 ± 0.2 ± 1.2 


iV° bs 


6 


3 




(1.5t™) x 10' 


5 (o.i4±S:S±£:S) x io~ 5 


Limits 


(> 0.4, < 3.7) x 10" 


" 5 < 8.1 x 10" 5 


B[B -> KvV) 


(i.4±5:S±S:l) x 10- 5 


Limits 


(> 0.2, < 3.2) x 10" 5 



90% CL of 

B(B+ -> K* + vv) < 6.4 x 10~ 5 , 
B(B Q -> K*°vV) < 12 x 10~ 5 , and (5) 
B(B -> K*vv) < 7.6 x 10~ 5 . 

The combined central value is B(B — > K*vV) — 
(3-8±l;g) x ID 

Since certain new-physics models suggest that en- 
hancements are possible at high sb values, we also 
report model-independent partial branching fractions 
(ABi) over the full sb spectrum by removing the low-ss 
requirement. The ABi values are calculated in intervals 
of s B = 0.1, using Eq. © (with the N° hs , iVf cak , iVf omb , 
and e^ lg values found within the given interval) multiplied 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 




FIG. 5: The central values (points with la error bars) 
of the partial branching fractions ABi versus sb, for (a) 
B+ ->■ K + vv, (b) B° -> K°vV, (c) B+ -> K*+vv, and (d) 
The subplots show the distribution of the fi- 
nal signal efficiencies within each sb interval (histogram with 
error bars), and the value of the total signal efficiency after 
applying all selection criteria without any requirement on sj 
(dotted line). The partial branching fractions are provided 
only within the intervals that are unaffected by the kinematic 
limit at large sb- 




FIG. 6: The constraints at 90% CL on e and r\ for sensitivity 
to new physics with right-handed currents. The B — > KvV 
(diagonal shading) and B — > K*vV (grey shading) excluded 
areas are determined from the upper and lower limits of this 
B — > K^*'vV analysis (solid curves), and from the most- 
strin gent upper limits from previous semileptonic-tag anal- 
yses [13, [l8|] (dashed curves). The dot shows the expected 
SM value. 



handed currents from new physics, such as non-SM Z° 
penguin couplings, could produce non-zero values. Using 
the parameterization from Ref . [l| , 



\ci m \ ' ^~|C£I 2 + W u 

the B — > K*vv upper limits from this search improve 
the constraints from previous searches on the Wilson- 
coefficient parameter space, as shown in Fig. [6] The B — > 
KvV lower limit provides the first upper bound on 77 and 
lower bound on e. These constraints are consistent with 
the expected SM values of e — 1 and 77 = 0. 



by the fraction of the signal efficiency distribution inside 
that interval. Figure [5] shows the partial branching frac- 
tions. The signal efficiency distributions are relatively 
independent of sb, which are also illustrated in Fig. [5] 
To compute model-specific values from these results, one 
can sum the central values within the model's dominant 
interval(s) (with uncertainties added in quadrature) and 
divide the sum by the fraction of the model's distribu- 
tion that is expected to lie within the same sb intervals. 
These partial branching fractions provide branching frac- 
tion upper limits for several new-physics scenarios at the 
level of lfr 5 . 

The B — > K^vv decays are also sensitive to the short- 
distance Wilson coefficients \C V L R \ for the left- and right- 
handed weak currents, respectively. These couple two 
quarks to two neutrinos via an effective field theory point 
interaction j^- Although \C R \ = within the SM, right- 



RESULTS FOR cc -»• vV 

In the search for cc — > i/U, Fig. [7] shows the m V v 
distribution of the observed data yields, expected back- 
ground contributions, and SM signal distributions. Ta- 
bles I VII and IVIII summarize the background contribution 
values and signal efficiencies within the J/ip and ^p{2S) 
invariant mass regions. The Tables also report the com- 
bined branching fraction central values and the branch- 
ing fraction upper limits at 90% CL for J/t/j — > vv and 
ip(2S) —> vv. The signal efficiencies account for the 
B -> J/ip and B -> K^ip(2S) branching fractions 
and their errors, which are taken from Ref. 29]. The 
data yield is consistent with zero observed cc — > vv sig- 
nal events in all channels. 

The combined upper limits for the charmonium 
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TABLE VI: Expected J/ip — > vv background yields Nf eak and N bks , signal efficiencies e^ lg , number of observed data events 
N° hB , and the resulting branching fraction central value and upper limit at 90% CL, all within the m„- invariant mass region 
corresponding to the J/tjj mass. Uncertainties are statistical and systematic, respectively. The Nf olnb yields are calculable as 



Channel 


K+ 


K° 


K *+ K+n° 


— > vv 

K*+ K° s n+ 


K*° -»• K + 7T" 




jypcak 


0.4 ±0.2 ±0.0 


0.7 ±0.3 ±0.1 


0.8 ±0.3 ±0.1 


0.4 ±0.2 ±0.0 


2.6 ±0.5 ±0.3 


0.6 ±0.2 ±0.1 


jybkg 


0.5 ±0.2 ±0.0 


0.7 ±0.3 ±0.1 


0.8 ±0.3 ±0.1 


0.8 ±0.3 ±0.0 


2.8 ±0.5 ±0.3 


0.6 ±0.2 ±0.1 


ef« (xlO- 8 ) 


95.3 ±4.4 ±5.5 


19.3 ±1.0 ±2.1 


20.9 ± 1.5 ± 1.7 


12.4 ± 0.8 ± 1.0 


36.2 ± 1.9 ±4.0 


1.8 ±0.2 ±0.2 


N° hs 


1 





1 








1 


B(J/i> ->• vv) 
Limit 






(0.2± 2 l±° i) x 10" 3 
< 3.9 x 10" 3 







TABLE VII: Expected tp(2S) — > w background yields jVP eak and N t s , signal efficiencies e^ lg , number of observed data events 
N° , and the resulting branching fraction central value and upper limit at 90% CL, all within the rrw invariant mass region 
corresponding to the tp(2S) mass. Uncertainties are statistical and systematic, respectively. The iV, comb yields are calculable 

Arbkg jvrpcak 



Channel 


K+ 


K° 


V(2S) - 
K *+ ^ k+tt 


> VV 

K*+ -> K° s ir + 


K*° -> K+tt- 




jypcak 


1.4 ±0.4 ±0.1 


0.6 ±0.3 ±0.1 


1.4 ±0.4 ±0.1 


1.0 ±0.3 ±0.1 


3.5 ±0.7 ±0.3 


0.6 ±0.2 ±0.1 


N hks 


1.6 ±0.4 ±0.1 


0.7 ±0.3 ±0.1 


1.4 ±0.4 ±0.1 


1.5 ±0.4 ±0.1 


3.9 ±0.7 ±0.3 


0.6 ±0.2 ±0.1 


ef e (xlO" 8 ) 


57.2 ±3.5 ±3.3 


13.1 ± 1.2 ± 1.4 


8.1 ±1.7 ±0.7 


4.9 ± 1.1 ±0.4 


14.2 ± 1.2 ± 1.6 


0.6 ±0.1 ±0.1 


iV° bs 


3 


1 


1 


3 


5 


1 


B(ip(2S) -> vv) 
Limit 






f5 6 +7 - 4+1 ' 6 

W- D -4.6-1.4 
< 15.5 X 


) x 10" 3 
10~ 3 







branching fraction values are determined to be 



B{J/ip -± vv) 
B(J/ip ->• e+e") 

B(ip(2S) -» vv) 
B(ip(2S) -> e+e-) 



< 6.6 x 10" 2 and 



< 2.0, 



(7) 



where B(J/jl> — » e + e~) and B(ip(2S) —> e + e~) are taken 
from Ref. [29[. With the addition of a new-physics J7 
boson, these ratios would be proportional to |/ c yc x>¥ ,|, 
where c XjV and / c y are the U couplings to the LDM 
particles \ or ip and to the c-quark respectively fl3j |. 
The J/ip decay ratio yields upper limits at 90% CL of 
\fcvCx,v>\ < (3.0,2.1,1.5) x 10~ 2 for spin-0, Majorana, 
and Dirac LDM particles respectively. These limits are 
comparable with those obtained by BES for J/ip — > vv, 
in ip(2S) -> tt+tt-J/^ 



a data sample of 471 x 10 6 BB pairs. We observe no sig- 
nificant signal in any of the channels and obtain upper 
limits at the 90% CL of B(B -> Kvv) < 3.2 x 10~ 5 , 
B(B -> K*vV) < 7.9 x 10~ 5 , B{J/ip -> vv) < 3.9 x 10~ 3 , 
and B(ip(2S) -> i/F) < 15.5 x 10~ 3 . The branching frac- 
tion central values and upper limits are consistent with 
SM predictions. We report B — > K^*'uv branching frac- 
tion limits in Tables IIVI and [Vj and cc — > vv branch- 
ing fraction limits in Tables IVII and IVII1 These results 
include the first lower limit in the B + — > K + vv de- 
cay channel, the most stringent published upper limits 
using the hadronic-tag reconstruction technique in the 
B° -> K°vv, B+ -> K*+vv, and B° -> K*°vv chan- 
nels, and the first upper limit for tp(2S) — > vv. We also 
present partial branching fraction values for B — > KS*' vV 
over the full sb spectrum in Fig. [5] in order to enable ad- 
ditional tests of new-physics models. 



SUMMARY 

In conclusion, we have searched for the decays B — » 
Kvv and B — > K*vv, as well as J/ip — > i/p and 
ip{2S) W via B -> K^J/ip and S ->■ K^ip(2S), re- 
coiling from a hadronically reconstructed i? meson within 
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FIG. 7: The m„- = sj SBm 2 B distribution for (from top to 
bottom) B + -¥ K + cc, B° -> K°cc, B+ -> if*+cc, and 
B° — > K*°cc events after applying the full signal selec- 
tion. The expected combinatorial (shaded) plus mES-peaking 
(solid) background contributions are overlaid on the data 
(points). The signal MC distributions (dashed) are normal- 
ized to B(cc — > vT>) values of 2% for the K + channel, 10% for 
the K° channel, and 5% for the K* channels. 
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